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SummaRT 


'because  or  excessive  computation  time,  solving  the 
parabolic  equation  in  (iiqner  dimensions  by  means  or 
implicit  Finite  difference  schemes  seems  to  Oe 
u.ipract teal  even  it  me  scneme  is  unconditionally 
jtaOle.  To  economise  tne  computation  time  and  computer 
storage,  a  staple  explicit  rimte  difference  scneme  is 
intruuuceu  for  the  solution  of  the  parabolic  equation  of 
the  ‘Penrod i nyer  type.  Tms  explicit  scneme  involes  five 
spatial  points  and  is  conditionally  stable  by  V) 
introducing  an  additional  dissipacive  term.  The 
complete  tneory  with  respect  to  the  stability  is 
proved.  An  application  to  a  tnree -dimensional  ocean 
acoustic  propagation  problem  is  included  to  demonstrate 
it,  val id! ty. 
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For  a  more  general  expression,  «e  can  include  the  low 
oroer  terms  to  give 
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As  an  application,  a  oneway  ocean  acoustic  souno 
propagation  In  three  dimensions  is  represented  oy 
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an/  ufiysicai  prooltf/i*  result  in  tne  real 
jgpiicatiun  ur  parahol il  equal i uns .  A  familiar 
re|jresentat we  parabolic  equation  is  the  neat  equation 
.jitn  real  coefficients*  *  nuniher  or  applications  (other 
tnji.  i»eat  cuuuudionj  iri^e  m  tne  area  of  quantum 
:ieci. jnics,  plaui.i  pnyika,  optics,  seismology,  ocean 
acoustics*  etc.  .  ij,  juu  result  in  a  form  of  parabolic 
equation  m  in  cixuiplex  coefficient*.  A  familiar 
representative  parauul ic  equation  with  complex 
coefficient*  is  the  >cnruumqer  equation.  Fur 
aiscussivin,  the  theory  or  a  new  stable  explicit  finite 
oirrerence  scneme  is  fell  as  i  real  application  are 
s.nuten  to  oeai  *itn  the  >cnr«jOinger  equation  of 
kj  1 1 1  -q».  nens  ions  mi  tne  form 
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where  k0  is  4  reference  wavemunDer  ana  n(r,*,jr)  is  tne 
three-dimensional  index  of  refraction,  wnicn  is  defined 
as  a  ratio  of  a  reference  sound  speed  to  a 
three-dimensional  sound  speed.  Eq.  (3)  is  in 
three-dimensional  cylindrical  coordinates  [2]. 


A  solution  exists  [3]  fur  Eq.  (3)  that  uses  an 
unconditonaf iy  stable  implicit  finite  dirrerence  scheme, 
winch  discretizes  Eq.  (3)  Oy  means  or  central  finite 
differences  for  Both  i  ana  »  derivatives.  Then  tne 
Cranx-Ntcolson  scneme  is  applied  to  formulate  a  large 
system  of  spared  matrix.  Tins  system  was  solved  Dy  a 
Tale  University  [3]  preconditioning  sparse  technique. 
Results,  proouce u  Oy  tne  Crank -Nicol son  scneme  [2J  are 
reasonably  accurate,  however,  oue  to  the  step-by-steu 


iteration  to  solve  the  system,  excessive  computer  tune 
was  required.  This  motivated  us  to  develop  a  more 
economical,  stable  explicit  finite  difference  scheme. 

In  the  sections  to  follOM,  tne  main  discussion  is  on  tne 
introduction  of  a  conditionally  stable  explicit  finite 
difference  thoroughly  examining  its  consistency, 
stability,  and  convergence.  A  theorem  to  describe  tne 
stability  of  this  nee  scheme  is  developed  and  proved, 
following  the  theoretical  section,  we  use  a 
three-dimensional  acoustic  wave  equation  arising  from 
the  application  of  underwater  wave  propagations  as  a 
test  case  to  examine  the  validity  of  the  theory.  We 
examine  the  accuracy  and  speed  of  tne  tneory  by 
comparing  it  with  the  solution  produced  by  the 
Crank-hicolson  scheme.  As  a  physical  illustration  of 
the  three-dimensional  problem,  a  plot  is  Included  to 
describe  intensity  effects  of  the  three-dimensional 
ocean  wave  propagation. 


A  STABLE  EXPLICIT  SCHEME  FOR  HIGH  DIMENSIONS 
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In  Eq.  (8),  j  represents  a  multi-index  (Jj,  . 

J  ).  D?  .  is  tne  second-order  centered  difference 

m  J.4 

operator  with  respect  to  j,  and  n^is  the  corresponding 
mesh  size. 


THEOREM:  If  scheme  (8)  is  used  to  solve  Eq.  (7),  the 
™  scheme  is  stable  if  and  only  if  •  <  0  and 


Chan,  Shen,  and  Lee  [4]  discussed  the  solution  to  a 
model  Schrodinger  equation,  i.e. , 


The  least  restrictive  stability  constraint  is 
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by  tne  finite  difference  scheme 
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and  is  obtained  when  a  «  -1/4  and  a  -  1/4. 

Tne  proof  appears  in  Its  entirety  in  reference  1 
and  is  outlined  below. 


which  IS  UNSTAULE  where  X  .  ar,  h  •  az. 


PROOF:  For  economy  In  writing,  define 


As  a  consequence,  a  number  of  stable  explicit 
schemes  were  introduced  [4]  to  solve  tne  parabolic 
equation  of  the  Scnrbdinger  type.  In  this  paper,  a 
scheme  is  selected  for  application  and  replaces  scheme 
(5)  by  introducing  «  dissipative  tarm,  which  is  added  to 
scheme  (5)  to  give 
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wnere  a  and  a  are  determined  to  be  a  •  -1/4,  »  •  1/4  for 
least  resitnctive  stability  condition.  As  a 
generalization  of  tne  scheme  lb)  to  the  Schrodinger 
equation  of  nign  oroer  [1],  consider  the 
multi -dimens lonal  Schrodinger  equation 
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The  amplification  factor  R  can  ba  determined  to  be 

m  , 

R  ’  1  6Jri  "(  la  *  '•)<i|). 
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The  stability  requires  that  ||w||  <_  1.  After  some 
simplification,  the  stability  condition  can  be  written  as 
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»nere  tne  o‘s  are  assumed  to  nave  tne  same  sign, 
without  loss  of  generality,  we  assume  o  >  u  tor  «  1, 
2,  ....  m.  we  consider  the  natural  extension  of  scheme 
(61  taxes  the  form 


Let  p  w  (nj,  hj,  ....  and  define 
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»«  consider  twu  cases:  a  _>  1/4  and  g  <  1/4. 

CASE.  1:  *  >  1/4. 
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lb:  In  02, 
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From  Iqs.  (Ill)  and  (111,  It  can  be  verified  that  for  g  _> 
1/4,  the  stability  condition  is 


x  <  min  (-  a/|8[a2  *  l»  *  -Jl2  3  ^  Fj) 


(12) 


CASt  11:  g  <  1/4. 

Clearly  Oj  is  empty  and  0  ■  02-  It  is  seen  that 
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Tne  general  stability  condition  is  therefore  (12). 
Clearly,  we  must  nave  a  <  II.  To  choose  a  and  g  such 
that  the  stability  condition  is  trie  least  restrictive, 
we  must  take  «  -  1/4  so  that 


k  <  min 
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To  maximize  tne  right -hand  side  above,  we  take  a  *  -1/4, 
wmch  gives 

k  — — — g —  ,  establishing  tne  result. 
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In  tne  tnree -dimensional  ocean,  a  class  of  sound 
wave  propagation  problems  can  De  represented  by  a 
parabolic  equation  of  tne  Penrod  I  riyer  type  [2J.  For 
prescribed  environmental  conditions,  an  application  of  a 
three-dimensional  problem  in  sector  can  be  shown  as  in 
Figure  1  for  its  sector  region  of  progagation,  wnere 
r0  r  r,„,  i)‘  >  •  (  l’,  and  d  <_  t  <.  Hldm.  In 
actual  TimuUlion  the  sector  is  taken  to  be  -20 ’  <  «  < 
20*. 


t.'. 


Figure  1:  Sector  Region  of  Propagation 

An  exact  solution  u(r,»,z)  has  been  obtained  [2] 
and  takes  the  fora 

1  XT  <10) 

u(r.e.i)  -  sin(Q  z)  e’“*  e  0  , 

which  satisfies  the  three-dimensional  parabolic  equation 

(3). 

The  computation  speed  among  explicit  finite  schemes 
[1]  and  an  implicit  finite  scheme  [2]  using 
two-dimensional  as  well  as  three-dimensional  examples 
has  been  examined  by  Chan  et  al.  [5],  The  same 
three-dimensional  problem  with  known  exact  solution  was 
used  by  Chan  et  al .  [5]  to  examine.  In  particular,  the 
computation  speed  between  each  explicit  scheme,  as 
aescribea  In  [1],  against  the  implicit  finite  difference 
scneme,  as  described  in  [2].  Their  findings  snow  a  more 
favorable  compuatlon  speed  for  the  explicit  scheme  tnan 
the  implicit  scheme.  Me  extend  their  study  to  some 
three-dimensional  effects  using  the  explicit  scneme, 
expressed  by  tq.  (6). 

In  the  application,  the  Q  is  assigned  to  be  »/U)0 
and  the  modal  Index  m  Is  taken  to  De  3.  The  source  is 
placed  at  SUm  below  the  surface  and  propagates  the  sound 
in  a  regular  three-dimensional  cylindrical  region.  The 
propagation  is  required  to  reach  the  maximum  range  at 
SbUn  where  we  can  see  three-dimensional  effects.  Me 
limit  tne  propagation  to  a  sector  of  40*  (l.e.,  from 
-20’  to  ♦20*)  and  centered  at  the  origin  (0,0,0).  For 
simplicity,  the  three-dimensional  sound  speed  c(r,»,z) 
is  taken  as  a  constant  and  the  medium  is  assumed 
homogeneous.  Initial  boundary  values  are  generated  frexn 
tne  exact  solution. 

Since  our  numerical  results  produce  field  intensity 
information  at  all  receiver  depths,  we  can  output 
contour  plots  for  each  angle  ».  Figure  2  presents  a 
contour  plot  of  energy  flow  at  *  •  0*. 
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Figure  2:  Contour  plot  of  Field  Intensity 

The  plot  was  produced  by  the  existing  three-dimensional 
Crank -hi  col  son  scheme  in  conjunction  with  a  Tale  sparse 
tecnnique.  The  accuracy  of  the  results  have  been 
discussed  in  [2].  The  same  calculation  was  performed  by 
the  explicit  scheme  ( Eq.  ( 8) )  using  the  same  range  step 
size  (O.OOlm)  as  used  by  the  Crank -Nicolson  scneme.  The 
explicit  scheme  solution  produced  results  very  close  to 
Crank -Nicolsun's,  thus,  generated  the  same  plot  curves 
as  described  in  Figure  2.  however ,  the  advantage  of  the 
explicit  scheme  is  tne  CPU  time  required  for  the 
complete  computation,  which  is  approximately  3.S  times 
faster  tnan  the  Crank -Nicoloson  scheme  for  achieving  tne 
same  accuracy. 
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CONCLUSION 


we  nave  introduced  an  explicit  finite  difference 
scneme  to  solve  the  Scnrodlnger  equation.  This  scheme 
was  developed  to  De  conditionally  stable.  Numerical 
results  demonstrated  Us  accuracy  and  agreement  not  only 
witn  tne  Crana -Nicolsun  scheme  out  also  with  the  known 
exact  solution.  It  is  expected  that  If  this  scheme  Is 
implemented  in  a  vectorized  computer,  its  storage, 
implementation,  and  computation  time  advantages  would 
becume  evident.  We  snowed  only  one  of  the  explicit 
scnemes  we  nave  developed,  we  believe  other  explicit 
scnemes  we  have  introduced  [4]  may  nave  equal  advantages 
over  implicit  schemes. 


